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Abstract 

Few years ago, application of the mean field Bethe scheme on a given system 
was shown to produce a systematic change of the system intrinsic symmetry. For 
instance, once applied on a ferromagnet, individual spins are no more equivalent. 
Accordingly a new loopwise mean field theory was designed to both go beyond 
the one site Weiss approach and yet preserve the initial Hamitonian symmetry. 
This loopwise scheme is applied here to solve the Triangular Antiferromagnetic 
Ising model. It is found to yield Wannier's exact result of no ordering at non-zero 
temperature. No adjustable parameter is used. Simultaneously a non-zero critical 
temperature is obtained for the Triangular Ising Ferromagnet. This simple mean 
field scheme opens a new way to tackle random systems. 
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1 Introduction 



Collective phenomenon are rather difficult to solve exactly. Up to date, only some one 
dimensional problems and the square zero field Ising model allow an exact analytical so- 
lution [1]. To compensate this situation, a rich family of approximate methods has been 
developed over the last one hundred years. The most powerful one being the renormal- 
ization group techniques [2]. 

At start was the Mean Field Theory (MFT). It offers a very practical and simple 
tool to solve most collective phenomena [1]. While it is completely universal and generic, 
associated quantitative results are unusually poor. In particular critical temperatures 
and exponents are rather far from exact estimates [2]. Sometimes even the order of the 
transition may be wrong like for the instance in the Potts model [3]. 

The crudest and most simple version of MFT is the 1907 Weiss pioneer model [4]. 
It reduces the infinite number of fluctuating degrees of freedom down to one, Si, which 
couples to homogeneous mean field degrees of freedom m. The thermodynamics is then 
solved calculating the associated partition function from which the self-consistent equation 

< Si >= m (where < ... > means thermal average) is derived. 

In the case of Ising systems with q nearest neighbor interactions, Weiss theory gives 

< Si >= tanh(Kqm) where K = (3J, J is the exchange coupling, (3 = -j^f, &b is the 
Boltzman constant and T is the temperature. Associated critical temperature is K c = ^. 
At odd with the known exact result a phase transition is obtained at d — 1 (q — 2) [1]. 

From there it took 28 years before Bethe improved the Weiss model [5]. Instead of 
just one fluctuating spin, he considers a cluster of fluctuating spins with a central one 
and its nearest neighbors. The main achievement of the Bethe approximation is to yield 
the exact result of no ordering at one dimension. However, critical temperatures given by 
K c = tanh _1 (4r), 

are not much better than from Weiss model. Critical exponents stay 
unchanged. Latter on, using computer capabilities, larger size fluctuating clusters have 
been considered to obtain better critical temperatures [6]. 

However, few years ago the Bethe cluster scheme was showed to systematically change 
the system intrinsic symmetry [7]. Starting from a system with equivalent sites like for 
instance a square Ising Ferromagnet, It ends up making individual sites unequivalent. At 
this stage it is worth to stress that an approximation can be very crude and yet not wrong 
as long as it preserves the intrinsic symmetry of the problem. Otherwise it does change 
the physics of it. We are not talking here about a symmetry breaking of the higher phase 
symmetry as it occurs in a usual phase transition but of a change of the symmetry of the 
disorder phase itself. 

On this basis the challenge was to find out if it is indeed possible to build a MFT which 
considers more than two fluctuating spins, yet preserving the initial lattice symmetry. 
Indeed, Galam showed it is possible using a loopwise scheme (LWS) which articulates on 
finite-size one-dimensional closed loops [7] . Paving the whole lattice with these loops, half 
of them are kept fluctuating while the other half is averaged out with mean field degrees 
of freedom. The scheme is illustrated in Figure 1 for the square lattice. 

The LWS is a generic model. It was applied to a large class of ferromagnetic systems 
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Figure 1: The loopwise scheme in the square case: si, s2, s3, s4 are the fluctuating spins 
while ml, m2, m3, m4 are mean field averages 
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on Bravais lattices [7, 8]. It reproduces the exact result of no ordering at one dimension. 
Moreover, for Ising hypercubes, it exhibits a lower critical dimension d\ for long range 
ordering which is equal to the Golden number d t = 1+ 2 v/ ^ . However critical exponents are 
unchanged from Weiss model. 

On this basic, to determine the range of validity of this new LWS, it is of interest 
to check if it can yield new properties which are out of reach of previous mean field 
theories like frustration. For instance, when applied to the fully-frustrated Triangular 
Ising Antiferromagnet (TIA) most MFT predict a transition at a non-zero temperature 
while an earlier exact argument by Wannier proved no symmetry breaking occurs at any 
non-zero temperature [9]. 

Few years ago, to bridge this difficulty Netz and Berker introduced the hard spin 
recipe [10]. It combines a mean field calculation with some Monte Carlo sampling. When 
applied to the TIA, it yields the correct result of no ordering at T ^ 0. Later Banavar et 
al suggested that the Monte Carlo sampling could be reproduced by expanding all possible 
products of the 6 nearest neighbors spins of the "exact spin" but it was then disproved 
by Netz and Berker[ll]. 

More recently focusing on the TIA, Monroe approximated the triangular lattice with 
a Husimi tree built up of triangles [12]. It then allows to include properly frustration to 
get a correct phase diagram. However an Huzimi tree is not a triangular lattice. 

In this paper we apply the very simple LWS to the fully frustrated Triangular Ising 
Antiferromagnet (TIA). The Wannier exact result is recovered [10] and a transition is 
found at T = 0. The following of the paper is organized as follows. Section 2 deals with 
the frustration effect. In Section 3 the LWS is presented. The TIA is solved analytically 
in Section 4 using the LWS. In Section 5 using the same Equations, the Triangular Ising 
Ferromagnet (TIF) is also solved. Some possible applications are mentinned in the Last 
Section. 

2 The Frustration effect 

Frustration is a major ingredient of many physical systems. It results from the impos- 
sibility to minimize simultaneously all pair interactions. In turn it makes the ground 
state highly degenerate [9]. Frustration effects may arise from either quenched disorder 
or topological constraints. 

Random bond spin glasses are the archetype of frustration produced by disorder. The 
random distribution of quenched competing interactions generates analytical difficulties 
in calculating the thermodynamic functions. In particular to average the disorder over 
the logarithm of the partition function is yet a real theoretical challenge. Usual mean field 
treatments failed to incorporate simultaneously frustration and quenched randomness. 

On this basis the TIA has the advantage of being fully frustrated without any disorder 
making the study of frustration itself more easy. It is therefore the perfect candidate to 
check the ability of a new scheme to deal with frustration. In addition an earlier exact 
argument by Wannier [10] has proved the absence of symmetry breaking at any non- 
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zero temperature for this system. At constrast most mean field like approaches produce 
wrongly some non-zero critical temperature. Along this line, Netz and Berker recipe [10] 
with Banavar et al reformulation [11] stand at odd. 

3 The loopwise scheme (LWS) 

The LWS was introduced few years ago to overpass the symmetry inconsistency of the 
Bethe scheme, yet retaining its physical feature of including several fluctuating degrees of 
freedom [7]. 

To implement the LWS on any lattice requires to single out two identical interpene- 
trating sublattices. Each element being composed from a closed compact loop of degrees 
of freedom. The shape and number of these degrees of freedom are determined by the 
lattice topology. It is the smallest closed linear loop. For instance in the square case (Fig. 
1) it includes 4 spins while for the triangular lattice (Fig. 2) 3 spins are involved. One of 
the sublattice is fluctuating and the other one is mean field. 

Both sublattices are coupled via nearest neighbor interactions. The problem is thus 
mapped onto decoupled one-dimensional closed fluctuating chains in external fields. The 
fields originate from the coupling to the mean field loops. At this stage an exact analytical 
calculation can be performed whatever the chain size is. It is worth to note no adjustable 
parameter is used. 

The LWS is a generic model. It was applied to a large class of ferromagnetic systems 
[7, 8]. Being built on using closed linear loops it should be well adapted to embody 
frustration effects [9]. 

4 Solving the Triangular Ising Antiferromagnet 

We now apply the LWS to the fully frustrated TIA. We first partition the triangular 
lattice into two interpenetrated triangular sublattices A and B. Thermal fluctuations are 
then ignored on the B-sublattice while preserved within the A-sublattice. These triangles 
are closed loops with no center (see Figure 2). 

All nearest neighbor (nn) plaquettes of a A-plaquette are B plaquettes, and vice versa. 
Therefore, on a given plaquette each spin has two nn spins of the same species (within 
the same plaquette), and four nn spins of the other species (belong respectively to three 
different nn plaquettes). 

Above breaking of the initial lattice symmetry makes the partition function calculable 
by decoupling the fluctuating triangles. The A-sublattice degrees of freedom can thus 
be integrated out in the partition function. The initial lattice symmetry will be restored 
latter using the usual mean field self-consistent constraint (Eq. (|j) below). 
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Figure 2: The loopwise scheme in the triangular case: si, s2, s3 are the fluctuating spins 
while 1, 2, 3 represent mean field averages ml, m2, m3 




4.1 Setting the Equations 

Given a A plaquette, we label the 3 fluctuating spins Si, S 2 , S3. We then introduce 3 
magnetizations mi, m 2 , m3 for corresponding B plaquettes (Fig. 3). The Hamiltonian 
then writes, 

H = -J(SiS 2 + S 2 S 3 + S3S1) - 5 J (Si(m 2 + m 3 ) + S 2 (m 3 + mi) + S3 (mi + m 2 )) , (1) 

where 5 = 2 accounts for the coupling to the B mean field plaquettes. From Eq. ([I]) the 
partition function is 

Z= J2 exp{-/3H} ; (2) 

Si=±l 

where % = 1, 2, 3. The three thermal average of Si, S 2 , S 3 are given by, 

<^>=4 E SieM-PH} • (3) 

We can thus write the associated three self-consistent equations, 

< Si >= rrii . (4) 

4.2 Looking for minima 

Indeed we are looking for minina of the free-energy which results from the partition 
function Z. It is then worth to stress not all solutions of Eq. (|j) are minimum. A 
criterium to make Eq. (§) a derivative of a function is to require its cross derivatives with 
respect to the m,j to be equal, i. e., 

d d 

< S i>=^ZT< S i >> ( 5 ) 



9mj 9m 
for i, j = 1, 2, 3 

Writing, S = (< Si >, < S 2 >, < S3 >), anc? m = (mi, m 2 , ms), the problem is now to 
find a set S = {m E R 3 ; S(m) = m}, such that there exists a function F obeying to, 

(m - S(m)) = dF(m) = . (6) 

To solve it, we rewrite thermal averages < Sj > as, 

< Si >= — stf (31,82,33) , (7) 

A Si=±l 

where 

f{si, s 2 , s 3 ) = exp {K(sis 2 + S1S3 + s 2 s 3 ) + 5K(mi(s 2 + s 3 ) + m 2 (si + s 3 ) + m 3 (si + s 2 ))} 
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Let cr G S3 be a permutation. Considering cr(m) = (m CT (i) , m CT (2) , w CT (3) ) we have, 

<r(Z(m)) = Z(<x(m)), Z(-m) = -Z(m) . (9) 

Writing X = exp K and X; = exp^Xmj, < Si > are rational fractions in (x^, X) and we 
have, 



D 



(10) 



D = (1 + T 3 2 )X 4 + T 2 + T X T 3 , (11) 
<*> = l-2 X * (ri + r3 ^ ) + X4 , (12) 



where 



Ti = xi + x 2 + x 3 , 
T 2 = xix 2 + Xix 3 + x 2 x 3 , (13) 
T 3 = Xix 2 x 3 , 

are the elementary symetric functions. Note D > 0, X > 0, Xj > and |< 5$ >| < 1. 

Solving first the K = case, we get immediately < Si >= and the solution is 
rrii = 0. We can then proceed assuming K 7^ 0. 

4.3 The most general solution mi 7^ 1712 7^ m3 

We can now solve the equations, starting with the most general case mi 7^ m 2 7^ m 3 . Eq. 
@ is equivalent to, 



d < Si > dxj _d< Sj > dxi 



dxj drrij dxi dirti ' 



(14) 



that is, 
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(15) 



Suppose first, three different values of mj. It makes X\ 7^ x 2 7^ x 3 since X 7^ 0, which 
in turn, solving Eq. (|i5|) implies, 
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and, 

X* = — ^— . (17) 
T 3 2 - 2 v ' 

In conclusion 

D = (1 + T 2 )X 4 + T 2 + TiTg, = (18) 

which is impossible since D > 0. Therefore, we can conclude that out of the three m,, 
two must be equal. We then suppose mi = m 2 ^ m 3 . 

4.4 The solution exhibits the symmetry mi = ^ ms 

From above calculation we restrict the minima search to the subspace of solution mi = 
m 2 7^ m 3- It implies x% = x 2 and, 

/C ^ ^ J 2(1 +2^3) 

< Si > - < S 3 > = -2 — (xi-xsjj (19) 

so it makes, 

expjSKmx) - exp(SKm 3 ) ^ Q Q 
mi — m 3 

which in turn makes X < 0. Let us define P = X2X3 and iV = x\x 3l it gives, 

P 3 - N 2 

" P(N 2 P + 2N 2 - 2P - 1) ' 1 ; 

and, 

< Si > + < S 2 > + < S 3 >= ^j. (22) 

If P — 1 or N — 1 then mi = m 2 = m 3 = which is not possible since we assumed above 
mi = m 2 / "23- So (P — 1)(N — 1) ^ 0. On the other hand, as K < 0, we must have, 

I — N 

>0, (23) 



mi + m 2 + m-3 
which makes (P - 1)(JV - 1) < but Eq. (P) gives, 

4 1 P 3 - N 2 

X = P(l -P) + (2 + P)(iV 2 - 1) < °' ^ 

which is impossible. In conclusion all three rrii must be equal. On this basis we now 
assume mi = m 2 = m 3 = m. 
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4.5 The solution is fully symmetrical with mi = = = m 

We have now proved the minima belong to the solution subspace defined by the symmetry 
condition mi = m 2 = m 3 = m. On this basis, writing Y = exp(5Km), Eq. @ becomes, 

m _ Mm1 _ (r 4 -i)((F 8 + y 4 + i)^ + n 



so we have either m = 0, or both, 



m m 



> 0, (26) 



and K > 0. We also deduce that |/^r(m)| < 1. 

As fx(—m) = — /if(m) it is enough to solve the case m > 0. We thus obtain X > 
1,Y > 1, K > 0, or m = 0. Computing the derivative gives, 

ff(m) = 5 Y 4 (3 Y*X* + (Y 16 + 4 F 12 + 4 F 4 + 1) X 4 + 3 F 8 ) 

M j 1 (y 4 + i) 2 ((y 8 -y 4 + i)x 4 + 3F 4 ) 2 ' 1 ; 

Y 4 (Y 4 — 1 1 

f'( m ) = -325 2 K 2 * ^ ^ ?SJr (m), (28) 

KK ' (y 4 + i) 3 ((y 8 -y 4 + i)x 4 + 3F 4 ) 3 v 7 v ; 

where 

9K (m) = (y 4 + l) 6 + (2y 8 + 7F 4 + 2)(r 4 + l) 4 (X 4 -l) 

+ (y 16 + 7Y 12 + 21 Y 8 + 7F 4 + 1)(F 4 + 1) 2 (X 4 - l) 2 

+9y 8 (r 8 + y 4 + i)(x 4 -i) 3 >o. (29) 

Therefore when m > 0, in addition to < /^(m) < 1, we have f' K {m) > and 
fxi. 171 ) < 0- These properties allow to conclude that, 

1. If f' K (0) < 1, is the only fixed point of fx- 

2. If f' K (0) > 1, fx has exactly three fixed points, 0, b, —b where — 1 < b < 1 
Computing then, 

f'rtrt ^ 3exp(4if) + l 

/x (0) = 2^ exp(4g)+3 , (30) 
it appears to be an increasing function of K. It makes, 

f'M = 1, (31) 

to have a unique solution Z^o. Moreover, if K < Kq then /x(0) < 1 and if K > Kq then 
fk(0) > 1. 
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4.6 The actual minima 



Looking for minima of Fk where, 

dF K {m) 
dm 



= m-f K (m), (32) 



depending on the value of K, two cases appear quite naturally for K > K and K < K . 
It shows the Triangular Ising both Anti and Ferromagnets are solved simultaneously. 

4.6.1 First case: K > K 

In this case, f'Kijn) = m has 2 solutions m = and m 2 = a where a is a positive function 
of K. Having, 

F' K {m = 0) = 1 - f' K (m = 0) < , (33) 

m = is a maximum for Fk- In parallel, 

F K (m = y/E) = F K (m = -yfi) = 1 - f K ( mi ) > . (34) 

Therefore m = y/a and m = —\fa are minima of Fk- They correspond to the Triangular 
Ising Ferromagnet symmetry breaking at low temperatures where K is the associated 
critical temperature. 

4.6.2 Second case: K < K 

Then the unique solution of fxijn) = m is m = 0. There, 

F K (0) = 1 - f' K (0) > 0, (35) 
so it is a minimum for F K . This case embodies indeed two different physical situations. 

1. The first range of positive K, < K < K , corresponds to the disordered phase of 
above Triangular Ising Ferromagnet. 

2. At the same time, the range of negative K (K < 0) corresponds to the Triangular 
Ising Antiferro magnet. For this system the unique solution is always m = for 
the whole range of temperatures T > 0. We conclude situation proving thus no It 
means no ordering occurs for the TIA at any non zero temperature. The Wannier 
argument is thus recovered [10]. 

4.7 A transition at T = 

From the exact Wannier solution the Triangular Ising Antiferromagnet is known to exhibit 
a phase transition at T = to an ordered phase with broken symmetry among the three 
sublattices. Accordingly we now examine what our scheme yield in the case K — > oo. To 
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solve the equations it is more convenient to rewrite Z and < Si > in terms of T = tanh(K) 
and ti = tanh(<5 Krrij) . We first note |< Si >| < 1 since |sinh(x) | < cosh(x). Then, once 
the rrii are fixed within [—1, 1], the condition K — > oo make the ij to go to either one of 
the three values —1, 0, 1. 

Computing < Si > in terms of ti and K for each one of the 27 possible limit values of 
the U set, we find 7 solutions for the m, which are respectively, 

m . = 0, i = 1,2,3. (36) 

and 



mj = m.,- = — rrik = ±1. (37) 

To determine the actual minimum at T = we compute the associated values for free 
energy F = — /c^Tlog Z. The first solution mi = m-i = = yields, 

F=~ log((6 + 2exp(4ZO) exp(-tf)) > 1 (38) 
and for mi = irt2 = 1, m3 = —1 we get, 

F = log (2exp(— ZT) cos(25i^)(3 + exp(4ff))) — > 1-25, (39) 

i\ A — > — oo 

making the solution mi = m 2 = l,m 3 = —1 the minimum. However from Eqs. (38) 
and (39) the two free energies of the ordered/disordered phases are expected to become 
equal only at some non zero temperature, a little bit above zero temperature, that is 
quite close to a critical point. It is coherent to the known result of a phase transition for 
the Triangular Ising Antiferromagnet at T = in agreement with the previous improved 
mean field theory by Netz and Berker [10]. 



5 The Triangular Ising Ferromagnet 

Coming back to the TIF, we can go further and evaluate the value of the critical tem- 
perature Kq. At this stage it is worth to notice, all above results are independent of the 
value of 5 which accounts for the coupling to the mean field loops. 

Since 1 < ^- — r— - — < 3, when K > 0, from Eq. (1311) we obtain, 

exp(AK) + 3 

h 2 K " 2 b (40) 

In addition, in the limit of large 5, we get, 

rz 1 2 3 29 1 

K = 7^-^ + ^~WJ77 + O « • (41) 



25 5 2 A5 3 24 S 4 \5 5 
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To get a numerical estimate of the ferromagnetic critical temperature Kq requires to have 
the 5 value. 

From Eq. (JJ) a straightforward arithmetic leads to 5 = = 2 since 2 nn are treated 
exactly within the fluctuating loop out of the 6 triangular n n. Plugging then, 8 = 2 
into Eq. ( |5TD yields Kq = 0.1772. It is rather far from the exact numerical estimate 
K c = 0.2746 [14]. In comparison, a usual mean field gives Kq = g = 0.1667, while for 
Bethe it is K = tanh _1 (|) = 0.2027. 

6 Conclusion 

In conclusion, we have showed that the very simple and generic mean field Loopwise 
Scheme, proposed by Galam [7], is able to solve exactly the Triangular Ising Antiferro- 
magnet. Without any adjustable parameter it recovers the exact Wannier argument of 
no ordering at T ^ and a transition at T = [10]. From the same Equations the Tri- 
angular Ising Ferromagnet is also solved simultaneously. A phase transition is obtained 
into a ferromagnetic phase at a non-zero critical temperature 

Moreover, contrary to the Bethe scheme, it preserves the initial lattice symmetry, yet 
going beyond the one-site Weiss approach. It also yields no transition for Ising hypercubes 
at d = 1 with a lower critical dimension of di = 1+ 2 . 

The Loopwise Scheme should allow a new solving of a very large class of physical 
systems, in particular random systems with frustration. For future work we consider to 
apply it first to the Triangular Ising Antiferromagnet in a finite field and then on the 
stacked 3D version of it. Application to the Random Field Ising model should also be 
done. 
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